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Abstract. We show that in addition to well known Bessel, Hermite-Gauss, 
and Laguerre-Gauss beams of electromagnetic radiation, one may also construct 
exponential beams. These beams are characterized by a fall-off in the transverse 
direction described by an exponential function of p. Exponential beams, like 
Bessel beams, carry definite angular momentum and are periodic along the 
direction of propagation, but unlike Bessel beams they have a, finite energy per unit 
beam length. The analysis of these beams is greatly simplified by an extensive 
use of the Riemann-Silberstein vector and the Whittaker representation of the 
solutions of the Maxwell equations in terms of just one complex function. The 
connection between the Bessel beams and the exponential beams is made explicit 
by constructing the exponential beams as wave packets of Bessel beams. 
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i We dedicate this paper to the memory of Edwin Power who recognized in his book the power 
of Riemann-Silberstein vector although he has not been aware of its early history. 
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1. Introduction 

Directed beams are very common forms of electromagnetic radiation, especially in the 
optical range. Several mathematical representation of such beams were introduced in 
the past. Among various beams, the Bessel beams and the Laguerre-Gauss beams play 
a distinguished role because they are characterized by a definite value of the projection 
of the angular momentum on the direction of propagation. Beams carrying angular 
momentum were studied theoretically and experimentally (a collection of papers on 
the optical angular momentum has been recently published A presentation of 

new examples of such beams seems to be worthwhile. In this paper we produce a new 
family of beams — the exponential beams. The exponential beams are exact beam- like 
solutions of the Maxwell equations with a definite value of angular momentum and 
with an exponential fall-off in the transverse direction. They can be placed halfway 
between the Bessel beams and the Laguerre-Gauss beams. Similarly to the Bessel 
beams, exponential beams have a definite value of the wave vector in the direction of 
propagation and like the Laguerre-Gauss beams they are confined in the transverse 
direction. The exponential beams are not monochromatic, but their spectrum may 
be highly peaked. They can be conveniently represented as wave packets of Bessel 
beams. In our analysis of exponential beams we shall employ the tool developed in 
our recent work (jjj — the Whittaker representation of the Riemann-Silberstein (RS) 
vector — that enables one to derive any solution of the free Maxwell equations from a 
complex solution of the scalar wave equation. We begin with a brief summary of this 
representation. 

2. Whittaker representation of the RS vector 

A natural tool in the analysis of the solutions of the Maxwell equations, in both 
classical and quantum theories, is a complex vector F, 



that carries all the information about the electromagnetic field. We named this vector 
in 0] the Riemann-Silberstein (RS) vector because it made its first appearance in [S] 
and was more extensively analyzed in The RS vector has been occasionally used 
in the past [7JIB] in classical electrodynamics, but we believe that its power lies also in 
bridging a gap between the classical and the quantum theory of the electromagnetic 



Owing to the linearity and the homogeneity of the Maxwell equations, the 
coefficient in the definition JTJ is not important and we shall often disregard the 
question of the overall normalization of F and use the RS vector in the form 
F = E + \cB. Two pairs of Maxwell equations expressed in terms of F reduce 
to just one pair 



Following Whittaker [HI Elj , we shall express the solutions of these equations by one 
complex function %(r,t) (in fact, Whittaker used an equivalent representation of the 
solutions of real Maxwell equations in terms of two real functions) 




(1) 



field 0. 



d t F(r,t) 



icV x F(r,t), V-F(r,t) = 0. 



(2) 




(3) 
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where x(r, t) obeys the d'Alembert equation 

(ld t 2 -A) X (r,£) = 0. (4) 
c l 

In the cylindrical coordinates (p,4>,z), that are useful in the description of beams 
carrying angular momentum, the components of the RS vector are 

F p \ ( d p d z + (i/cp)d^d t \ 

f U (i/ P )d <p d z -(i/c)d p d t \xM,z,t). (5) 

F, ) V -d 2 p -(l/p)d p -(l/ P 2 )dl J 
3. Exponential beams 

We shall construct exponential beams as wave packets of Bessel beams. In this 
manner we can automatically obtain a spectral decomposition of exponential beams, 
since every Bessel beam is monochromatic. The Bessel beams are obtained from the 
Whittaker representation J2J) or 10 by choosing the scalar functions in the form [H] 

Xl q±m (p,^z,t) = e- u 't u * t -<"- m *)j rn (q x p), (6) 

where uj q = cq = cy 1 q\_ + q z and a = ±1 determines whether the wave will have the 
right-handed or the left-handed polarization. Since Bessel beams form a complete 
set of solutions of the Maxwell equations, every solution can be represented as a 
superposition of Bessel beams. Owing to the linear dependence of the RS vector on 
X, the construction of such superpositions can be realized by simply superposing the 
corresponding functions x- I n order to obtain a wave packet of Bessel beams with 
a fixed angular momentum, we must include in this superposition only the Bessel 
functions with a given value of m. The general wave packet of this kind (for a fixed 
value of a) is 

/>oc />oo 

Xm a(p, 0, z, t) = e l ° m * / dq ± / dq z h{q Zl q^)e-^ c ^+^ t -^ J m {q^p), (7) 
JO J -oc 

where h(q z ,q±) is some weight function. For beam-like solutions, the longitudinal 
component q z should not be spread out too much and also it must be much larger 
than the transverse component q±. 

There are very few weight functions that allow for an analytic evaluation of the 
integrals over q± and q z in l [7| ). In addition to those treated in our previous work 

there are several weight functions that lead to exponential beams. The simplest 
example of a beam with an exponential fall-off will be obtained by fixing the value 
of q z and taking a superposition of Bessel beams with a varying q±. We shall only 
write down the formulas for a = 1. The case a = — 1 can be obtained by the complex 
conjugation. The expression (JJJ becomes now 

/>oo 

Xm?2 (p, c/>, z, t) = / dq ± / l (^)c- ic V^i i J m ( (7± p). (8) 

Jo 

Obviously, the beam constructed in this way will not be monochromatic — the weight 
function h(q±) determines the spectral characteristic of the beam. In order to show 
explicitly the distribution of frequencies, we change the integration variable from q± 
to oj = CyJ ql + q'j_. In addition, for mathematical convenience, we shall replace the 
frequency to by a dimensionless parameter w — u/c\q z \ = yl + (q±/qz) 2 and we 
rewrite (JSJ in the form 

Xmq ApA,z,t) =e^ + *^ J™dwg{w)e~ {c ^ wt J m (\q,\ P jvfl - l) , (9) 
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Figure 1. The normalized spectral weights g n {w) plotted as functions of w for 
cq z r = 50 (the lowest peak), cq z r = 100, and cq z r = 200. 



where g(w) is the spectral function. There are several spectral functions for which 
this integral can be evaluated analytically. The simplest example is (other functions 
are given in | Appendix A| ) 

g(w) = (w 2 - l) m/ V c fel™. (10) 

The parameter r has the dimension of time and it determines the spread of frequencies 
Aoj = 1/t in the wave packet. In figure ^ we show the normalized spectral functions 

9n(w) = , . y p-j — j— r, / dwg n (w) = 1 11) 

T(m/2 + l)K m / 2+1 / 2 {c\q z \T) J l 

for the three values of |<z z |cr. In order to better show the details, we have chosen in 
figure ^ broadband beams (small r). Most beams used in realistic experiments are 
quasi monochromatic. For such beams, the characteristic dimensionless parameter 
|q 2 |cr is several orders of magnitude larger. The corresponding spectral functions 
would have the form of very high and very sharp needles. The integration over w in 
|J3J, for g{w) given by ljTU)l . can be performed due to the relation (cf. Eq. 6.645.2 of 
Ref. [Til) 



f3 v K u+1/2 (v/a 2 +/3 2 ) 



dx (J - iy/ 2 e~ ax j v (/V^T) = t/- (q2 + /3 v 2) . /2+1/2 J (12) 



and it leads to the following result for \ 



H'lQz 



(13) 



where 

s = VP 2 -c 2 (t-ir) 2 . (14) 
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The Macdonald functions K m+ 1/2 (a?) for integer m reduce to ordinary exponentials 
multiplied by polynomials in the inverse powers of the argument 

ts , \ ~x\-^ (m + k)\ 1 

^W(-) = V ^ g4^)T W - (15) 

Thus, the functions Xmq z describe beams with an exponential fall-off. Indeed, for large 
\q z \p, but at the fixed values of all the remaining variables, these functions decrease 
exponentially 

Xmq ,M,z,t) ~ ji™4+***)%JL!L. (i6) 

P^oo \q z \p 
The first three functions Xmq z (for m = 0, 1, 2) are 

Viz I s 

X2 „ (p, A - e«+«">^ (l + ^ + ^) ■ (170 
In the simplest case, when m = 0, the RS vector is 

F P \ ( -iq zP {s 2 + \q z \s 3 ) \ pifei.-l9.la 

= ic(t-ir)p(3 + 3| & | S + (| fe | S ) 2 ) 1 115 ' ( lg ) 

y v 2 ( s2 + \v*\ s3 ) - p 2 ( 3 + 3 i^i s + (i^i*) 2 ) / lqzls 

For nearly monochromatic exponential beams, i.e. for very large values of the 
parameter r, as compared with p/c and t, the functions Xmq z have the form 

P — \qz\cT 

XmgM (x,y,z,t) ~ (g + iy^e-lff-M+H.* (19) 

|g z |(cT) m+i 

clearly indicating a straight vortex line of strength m along the z axis. This is a 
general characteristic of all beams of radiation endowed with angular momentum 3 . 
For m — 0, we obtain the standard plane wave. 

Having at our disposal the expansion © of the exponential beam into the Bessel 
beams, we may calculate the spatial Fourier representation of the exponential beam. 
To this end, we use the expansion of the plane wave into Bessel beams (cf. Eq. 8.511.4 
of Ref. [HI) 



e -ik-r 



e - ikz z (-i) m e im ^-^ J m (k±p) (20) 

m— — 00 

in the formula for the Fourier transform of Xmq z 

Xmq z (k±,ip,k z ,t) = /d 3 re _lfe " r Xmfe(p,0,^i)- (21) 



Upon substituting Xmq z , as defined by , into the formula l|21(l , we can easily perform 
the integrations over ip and z leading to delta functions. The remaining integration 
over p produces also a delta function since 

3 S(k q ) 

pdpJ m (k±p)J m (q±p) = (22) 
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Figure 2. The energy density F* ■ F (in arbitrary units) shown as a function of 
the distance p from the vortex center and time t. 



and the final formula reads 

k m 

Xmo, (ku<P, *.,*) = (-ir(2*) a *(*. - g.) {kz{m+1 j k , ± + fc| e^e-V^Hg(^).(23) 

This result can also be obtained directly, without the use of the expansion |(?JJ|, by 
evaluating the Fourier transform of the original expression (| 1 3|) for Xmq z , with the 
help of the formula 6.596.7 of Ref. [TT] 



4. Time evolution of exponential beams 

Exponential beams are not stationary — they evolve in time. At a given time, our 
exponential beams look like diffraction-less beams — their energy density and the 
Poynting vector do not change along the direction of propagation. In this respect 
they resemble Bessel beams. However, they exhibit a specific time dependence. In the 
transverse directions these beams first shrink and then expand. In our formulas the 
origin of time is chosen at the moment when the beam is maximally squeezed. This 
behavior is illustrated in figure [5] where we plot the energy density F* ■ F calculated 
for the RS vector l|18|l as a function of p and t. In order to describe the time evolution 
in quantitative terms we shall use as a measure of the transverse size of the beam the 
mean square radius of the energy distribution 

(x 2 + y 2 )(t) = J jdxdy(x 2 + y 2 )\F\ 2 /£ tot . (24) 

We prove in |Appendix B| that the time dependence of this quantity is simply quadratic, 
so that the formula for the time dependence reads 

(x 2 + y 2 ){t) = (x 2 + y 2 )(t = 0) + 2t 2 £ z /£ toU (25) 
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Figure 3. Time dependence of the mean square radius for m = 1 plotted for 
four values of the parameter r = 5, 10, 20, 40. Time and r are measured in units 
of 1/cgz and the radius is measured in units of 1/<J 2 . The larger the value of t, 
the wider the parabola. 



where £ z is the energy carried by the z components of the electromagnetic field and 
£tot is the total energy both evaluated per unit length in the z direction, 



£ z = 



J Jdxdy\F z \ 2 , £ tot = J Jdxdy\F\ 2 . (26) 

Several plots exhibiting the dependence on the parameter r are shown in figure [31 
Note that the more the beam is squeezed at t = 0, the faster it expands. 



5. Conclusions 



We have constructed a new class of exact beam-like solutions of the Maxwell equations 
with an exponential decrease in the transverse direction. They are characterized by the 
wave vector q z and the projection of the angular momentum m on the beam direction. 
These beams spread out in time but for a fixed value of t the energy distribution does 
not change along the beam direction. The exponential beams are not monochromatic, 
but they can be made nearly monochromatic by taking a large value of the parameter 

T. 
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Appendix A. Further examples of exponential beams 

In addition to the simplest case described before, there are several other examples of 
beams with an exponential fall-off in the transverse direction. They are derived from 
the following formulas involving integrals of Bessel functions (cf. 6.645.1 and 6.646.1 

of El) 



oo 

dx 



x ~ lY'L r ( a . /Ii T\ _ exp(-Va^TF) ( P 



dx(x 2 - l)- 1/2 e- ax J v (/V^ 2 ~ lj 



Upon choosing ff(to) in the form 

/ — 1 \ m / 2 

^ = (^tt) e " c '"'™' {L2) 

we obtain the following expression for an exponential beam 

e i{m<f,+q z z) e -\ qi \s m e i(mui+g^) e -|<?*| S / s + c ( T + i ^\™ 
Xm fe (P, <?>, 2, *) = | j— ; - ( , ... V n = i i > (1-3) 

\q z \s (s - c(t + it)) \q z \s \ P J 

where s was defined in (|14|) . 

The formula Ijl.lbjl requires several additional steps before leading to exponential 
beams. First, we shall put v = 1 and express the Bessel functions ix/2 an d K 1 / 2 by 
elementary functions. This leads to 

r °°dx (x 2 - l)-V2 e -a(*-i) Jt (V^ - !) = i (l - e Q ~V^+^) . (1.4) 

By differentiating this equation with respect to a and /? we may generate a whole 
family of new formulas. In particular, we obtain the formula 

JO 

dx {x 2 - l) m/2 {x + l)-^-™* J m (pVx 2 - l) 



■ "' ^ e ~y/ct 2 +f3 2 



that corresponds to 

g(w) = {w 2 - l) m/2 {w + ij-ie- ''-!" , (1.6) 

and it leads to another exponential beam. General characteristics and the time 
evolution of all exponential beams are very similar. The overall behavior, as discussed 
in the main text, is determined by the exponential factor exp(— |<? z |s). 

Appendix B. Time dependence of the mean square radius 

We can find the time dependence of the transverse mean square radius by 
differentiating twice the expression 



JJdxdy( 3 



s'+irW (2.i) 
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with respect to time and using the Maxwell equations. This is done in three steps. 
First, we use the formulas 

<9 2 |F| 2 = d 2 F*-F + 2d t F*-d t F + F*-d 2 F 

= c 2 (AF*-F + 2(V x F*)-(V x F) + F*-AF) 

= c 2 (A(F*-F)-2d t d J (F*F J )). (2.2) 

Note that the derivatives with respect to <f> and z drop out in the last line, because 
of the phase factor Q l (iz z + m< t>) in F. Next, we multiply the last line by x 2 + y 2 and 
integrate over x and y, to obtain after integration by parts 

d 2 J J dxdy(x 2 +y 2 )\F\ 2 = 4 J J dxdy\F z \ 2 . (2.3) 

The last step is to prove that the last integral does not depend on time. To show this, 
we again use the Maxwell equations, integrate by parts, and add and subtract terms 
with derivatives with respect to z 

d t J J dxdy\F z \ 2 = icj J dxdy (F*(d x F y - 8 y F x ) - (d x F; - d y F*)F z ) 

= ic J J dxdy (F* ■ (V x F) + {F;d z F x - F*d z F y ) . (2.4) 



Since this expression is real, we may add its complex conjugate and take a half of the 
sum to obtain 

| J J dxdy (F*-(V x F) — (V x F*)-F)) = — ^ / J dxdy\7 ■ (F* x F) = 0. (2.5) 

As a byproduct of this calculation, we also obtain the time independence of <?tot 
because the first integral (apart from 1/2) is the time derivative of £tot- Putting it all 
together, we have 

J J dxdy(x 2 + y 2 )\F\ 2 = J J dxdy(x 2 + y 2 )\F\ 2 + 2c 2 t 2 J J dxdy\F z \ 2 . (2.6) 

Dividing both sides of this equation by ftot, and taking into account the time 
independence of £ tot , we obtain (|25() . 
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